YITP-01-53 

UTCCP-P-108 

RCNP-TM1017 



Numerical study of 0(a) improved Wilson quark action 

on anisotropic lattice 

Hideo Matsufuru 1 , Tetsuya Onogi 1 and Takashi Umeda 2 

o ■ , 

O . Yukawa Institute for Theoretical Physics, Kyoto University, 

<N ; Kyoto 606-8502, Japan 

^ ! 2 Center for Computational Physics, University of Tsukuba, 

Tsukuba 305-8577, Japan 

30 June 2001 

>: 

Abstract 

C**** , The 0(a) improved Wilson quark action on the anisotropic lattice is investigated. We 

carry out numerical simulations in the quenched approximation at three values of lattice 
spacing (a^ 1 = 1-2 GeV) with the anisotropy £ = a a /a T — 4, where a a and a T are the 
spatial and the temporal lattice spacings, respectively. The bare anisotropy jp in the 
quark field action is numerically tuned by the dispersion relation of mesons so that the 
renormalized fermionic anisotropy coincides with that of gauge field. This calibration of 
Oh, bare anisotropy is performed to the level of 1 % statistical accuracy in the quark mass region 

below the charm quark mass. The systematic uncertainty in the calibration is estimated by 
comparing the results from different types of dispersion relations, which results in 3 % on 
our coarsest lattice and tends to vanish in the continuum limit. In the chiral limit, there is 

•'"Hi . 

k> , an additional systematic uncertainty of 1 % from the chiral extrapolation. 

Taking the central value 7f = Jp from the result of the calibration, we compute the 
light hadron spectrum. Our hadron spectrum is consistent with the result by UKQCD 
Collaboration on the isotropic lattice. We also study the response of the hadron spectrum 
to the change of anisotropic parameter, j F — ► 7^ + 5jp- We find that the change of 7^ by 2 
% induces a change of 1 % in the spectrum for physical quark masses. Thus the systematic 
uncertainty on the anisotropic lattice, as well as the statistical one, is under control. 

1 Introduction 

The anisotropic lattice is drawing more attention as a useful technique of the lattice QCD 
simulation in various physics such as the spectroscopy of exotic states, finite temperature QCD 
and the heavy quark physics. However, the advantage of having a fine lattice spacing in the 
temporal direction is obtained at the sacrifice of manifest temporal-spatial axis-interchange 
symmetry. Therefore improper use of the anisotropic lattice could lead to an unphysical result 
due to the lack of Lorentz symmetry in the continuum limit. This can be a serious problem 
for those physics in which the precision of the results is crucial. 

One way to avoid this problem is to tune the anisotropy parameters of the action by impos- 
ing the conditions with which the Lorentz invariance is satisfied for some physical observables. 
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In the Wilson plaquette gauge action there is only one anisotropic parameter in the action 
]!]]. Wilson loops are used to obtain the relation between the anisotropic parameter and the 
physical ratio of the lattice spacings in temporal and spatial directions, £ = a a /a T , where a a 
and a T are the spatial and the temporal lattice spacings [^, ^ [|, ||. On the other hand, not 
much is known about the quark action. This is because the previous works with the quark 
action on the anisotropic lattice have been devoted to the charmonium systems || |7|, §| || • In 
order to apply the anisotropic lattice to the systems containing the light quarks, one has to 
study how one can tune the parameter of light quark action in practical simulations. 

In this paper we study the 0(a) improved Wilson action on the anisotropic lattice using the 
quenched lattices with three lattice spacings, at fixed renormalized anisotropy, £ = 4. These 
scales covers the range of the spatial lattice cutoff a^ 1 =1-2 GeV. We first tune the bare 
anisotropy in the quark action numerically so that the renormalized fermionic anisotropy is 
equal to that of the gauge field by imposing the relativistic dispersion relation of mesons. This 
calibration is performed to the level of 1 % statistical accuracy in whole quark mass region 
below the charm quark mass. The extrapolation of tuned bare anisotropy to the chiral limit 
is performed by fitting to presumable forms, and causes additional systematic uncertainty of 1 
% at the chiral limit. The systematic uncertainty in the calibration from the lattice artifact is 
estimated by comparing the results from different types of dispersion relations. It results in 3 
% on our coarsest lattice, and tends to vanish in the continuum limit. Using the result of the 
calibration we compute the light hadron spectrum on anisotropic lattices, and examine how 
the uncertainty in the calibration affect the spectrum at the parameters of physical interest. 
It is found that the uncertainty of 2 % in calibration induces systematic error of 1 % in the 
spectrum. Thus the systematic uncertainty on the anisotropic lattice is under control, as well 
as the statistical one. We also show that the anisotropic lattices produce consistent results with 
those on the isotropic lattices, and that the Lorentz invariance of the simple matrix element is 
satisfied within errors. 

This paper is organized as follows. The next section describes the 0(a) improved Wilson 
quark action, which has been discussed in ju], |(|. The calibration procedures are discussed in 
Section ||. In Section ||, we perform the calibration of the bare anisotropy in the quark action. 
The systematic uncertainty induced by this tuning is fully examined. In Section |5| we apply 
the anisotropic lattice to the light hadron spectroscopy. In the last part of this section, the 
systematic uncertainty due to the anisotropy is again investigated in terms of the effect on 
hadron spectrum. The last section is devoted to our conclusions. 

2 Quark action on the anisotropic lattice 

2.1 Quark field action 

We employ the 0(a) improved quark action on the anisotropic lattice. The form of action 
has been discussed in Ref. JlQ] , which is the same as the Fermilab action [11| but defined on 



the anisotropic lattice. In this section we summarize the result which will be necessary in the 
following calculations. 

The quark action is represented in the hopping parameter form as 

S F = Y / ^)K(x,yMy), (2.1) 
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where n a and k t are the spatial and the temporal hopping parameters, r is the Wilson param- 
eter and ce and eg are the clover coefficients. In principle for a given K a , the four parameters 
K a /K T , r, ce and cb should be tuned so that Lorentz symmetry holds up to discretization 
errors of 0(a 2 ). 

On the anisotropic lattice, the mean-filed values of the spatial link variable u a and the 
temporal one u T are different from each other. The tadpole-improvement (12] is achieved by 
rescaling the link variable as Ui(x) — > Ui{x)/u a and U±(x) — > Ui(x)/u T . This is equivalent 
to redefine the hopping parameters with the tadpole-improved ones (with tilde) through K a = 
R a /u a and k t = k t /u t . We define the anisotropy parameter jp as 

IF = t. (2.3) 

At the tadpole-improved tree-level, and for sufficiently small quark mass, the anisotropy jf 
coincides with the cutoff anisotropy £ = a a /a T . 

In this work, we set the coefficients of the spatial part of the Wilson term and the clover 
coefficients as the tadpole- improved tree- level values, namely, 
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7, CE = 2 ' CB = ~^ 



r = 7> c E = m> °b = 7^ ( 2 - 4 ) 



and perform a nonperturbative calibration only for 7^ with the meson dispersion relation. 
It is useful to define k 

- = — - 2(7F + 3r - 4) (2.5) 

so that the bare quark mass in the spatial lattice unit, moo-, is expressed as 

m 0(7 = \ (i - 8) , (2.6) 

which is analogous to the case of the isotropic lattice. 

2.2 Dispersion relation of free quark 

In this subsection, we examine the dispersion relation of the free quark on the anisotropic 
lattice. First the tree-level relation of bare anisotropy with £ is derived from the condition that 
the rest mass and the kinetic mass coincide. Then, we discuss how the dispersion relation is 
distorted at the edge of Brillouin zone due to our choice, r = l/£ Q. 

From the action (|2.2|), the free quark propagator satisfies the dispersion relation 



P 2 + ("*> + l^P 2 " 2 
2(l + m + i^p 2 



cosh£( P ) = 1 + ",; ;••'" '. Yr tL . ( 2 - 7 ) 



where Pi = — sinpi, pi = 2 sin(pi/2) and mo = w-o<t/7f is the bare quark mass in the temporal 
lattice unit. Setting p = 0, Eq. ( |2.7D gives the rest mass 



Mi = £(0) = ln(l + m ). (2. 



On the other hand, the kinetic mass is defined and obtained as 
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Generally the rest mass ( |2.&j ) and the kinetic mass ( |2.9D are different. One can tune the bare 
anisotropy parameter jp so that it gives the same values for the rest and the kinetic masses 
[11]. Putting the rest and the kinetic masses equal, the anisotropy parameter ^p is represented 
with mo and r as 
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For small mo, 7_f is expanded in itiq as 
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The mo dependence starts with the quadratic term for r = l/£, therefore the dependence on 
the quark mass is small for sufficiently small mo- For example, let us consider the case of 
a T = 4 GeV, which corresponds to our coarsest lattice in the simulation. The charm quark 
mass corresponds to m$ ~ 0.3 and at this value jp is different from £ only 3 %. Up to this 
quark mass region, one can expect that the difference of 77? from £ is also small in the numerical 
simulation. This is examined in Section ||. 

With our choice of Wilson parameter, r = l/£, the action (2.2) leads the smaller spatial 
Wilson term for the larger cutoff anisotropy £. The question is how the contribution of the 
doubler eliminated by the Wilson term becomes significant for practical value of £. In the 
following argument on this subject, we only treat the case jp almost equal to £, i.e. the region 
of m$ sufficiently smaller than unity. Figure |] shows the dispersion relation ( |2.7| ) for several 
values of mo in the case of £ = 4, which we use in the numerical simulation. 

Let us examine the practical case, a^ 1 ~ 1.0 GeV, which corresponds to the lowest spatial 
cutoff of our three lattices. For the light quark mass region, mo = 0.02-0.05 corresponds to 
80-200 MeV, and roughly covers the mass region which we use in the hadron spectroscopy 
in Section ^. E(p) — E(0) rapidly decrease at the edge of the Brillouin zone, and the height 
at z = a/n is around 400 MeV. For two quarks with momenta p = ±cl/tt, additional energy 
of doublers is ~ 800 MeV, and is expected to affect not severely on the spectrum and other 
observables. For higher lattice cutoff, the situation becomes better. Then we regard the 
doubler contribution is sufficiently small on the lattice we use in the simulation, mo = 0.3 
roughly corresponds to the charm quark mass with a^ 1 ~ 1 GeV. In the case of heavy-light 
hadrons, such as D mesons and A c baryon, the scale of momentum transfered inside hadrons is 
of the order of Aqcd, and the same argument for light quark case holds. On the other hand, 
for the heavy quarkonium, the typical energy and momentum exchanged inside the meson are 
in the order of mv 2 and mv respectively ]i~3|| . For the charmonium, v 2 ~ 0.3, then typical 
scale of the kinetic energy is around 500 MeV. It seems not sufficiently smaller than the two 
doublers' contribution, and hence one need to choose larger lattice cutoff in the calculation of 
heavy quarkonium. 
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Figure 1: Dispersion relation of the free quark for £ = 4. 

3 Calibration procedures 

On anisotropic lattice, one must tune the parameters so that the anisotropy of quark field, £p> 
equals to that of the gauge field, £g : 



6K/3,7g;«,7f) =£g(/?,7g;k,7f) =£ 



(3.1) 



Since £g and £f are in general functions of both of gauge parameters (/?, jc) and quark param- 
eters, (k, Jf), the simulation with dynamical quarks requires to tune these bare parameters 
simultaneously. In the quenched case, however, this tuning is rather easy to be performed, 
since £g can be determined independently of jf- After the determination of £, one can tune 
jF so that a certain observable satisfies the condition (|3.1| ). In this work, we use the relativistic 
dispersion relation of meson, 



E 2 (p) = m 2 + ^ + 0( P % 
Zf 



(3.2) 



as our main calibration procedure. The energy and the mass of meson, E and m, are in the 
temporal lattice unit while the momentum is in the spatial lattice unit. The £p appears to 
convert the momentum in the spatial lattice unit into the temporal lattice unit, and considered 
as the fermionic anisotropy defined through this relation. With the condition £p = £, this 
condition satisfies that the rest mass and the kinetic mass equal to each other. For finite 
lattice spacings, above dispersion relation only holds up to the 0((p 2 ) 2 ) correction term. In 
the continuum limit, this higher order term in a would vanish and the relativistic dispersion 
relation would be restored. 

In the numerical simulation, we fit E 2 to the form Eq. ([T^) and obtain the value of £p for 
each input value of bare anisotropy 7^. Then we linearly interpolate £p in terms of 7^ and 
find out jp, the value of 7^ with which £p = £ holds. 



In order to estimate the systematic errors we also use the dispersion relation which corre- 
sponds to the lattice Klein-Gordon action ||, 



cosh E(p) - cosh E(p = 0) = 1 + -p^—p 2 - (3.3) 



Thus the comparison of these two calibration conditions typically shows the size of the lattice 
discretization errors. 

Expanding this expression in a, £k@ is related to £p as 



^G = ^U-^ + 0(a 4 )j. (3.4) 

The same input ^f gives smaller value for £,kg than £p, and therefore the tuned bare anisotropy 
7^ results in larger value in the former case. 

4 Numerical Results of Calibration 

The goal of this section is to determine the tuned bare anisotropy of quark field, 7^, at each 
fixed quark mass in the region from strange to charm quark masses. The reason for this choice 
of the quark mass range is that the simulation is easier which saves the amount of work in the 
exhaustive study for calibrations. Fitting the result as a function of the quark mass, we obtain 
j F to the statistical accuracy of 1 % level for the whole quark mass region below the charm 
quark mass including the chiral limit. 

Then we estimate the systematic uncertainties of 7^ which are mainly due to 0(aa) and 
0(a?) lattice artifacts. We also investigate how these systematic errors as well as the statistical 
error affect the meson masses, in the region m s < m q < m c . The response of hadron masses 
with respect to 7^ in the light quark mass region, m q < m s , needs additional care, and is the 
subject of next section. At the end of this section, we summarize the result of calibration. 

4.1 Simulation Parameters for the Calibration 

In this work, we use three lattices with j3 = 5.75, 5.95 and 6.10 and the renormalized anisotropy 
£ = 4. The value of jg corresponding to the desired value of £ has been studied in detail by 
Klassen [0], and we can use his relation of 7c and £ which is obtained in one percent accuracy. 
The statistical uncertainties are, otherwise noted, estimated by the single elimination Jackknife 
method with appropriate binning. The configurations are separated by 2000 (1000) pseudo- 
heat-bath sweeps, after 20000 (10000) thermalization sweeps at /3=5.95 and 6.10 (5.75). The 
configurations are fixed to the Coulomb gauge, which is particularly useful for the smearing of 
hadron operators. 

The lattice cutoffs and the mean-field values of link variables are determined on the smaller 
lattices with half size in temporal extent for (3 = 5.75, 5.95, and otherwise with the same 
parameters, while at /? = 6.10 the lattice size is 16 3 x 64. To obtain the lattice cutoffs, the 
static quark potential is measured with standard procedure. We adopt the hadronic radius ro 



proposed by Sommer |14| to set the scale. Following the method in Ref. [14], we determine the 
force between static quark and antiquark, as a function of 77, the interquark distance improved 
with the lattice one-gluon exchange potential form. Then we fit the values of force, containing 
the off-axis data, to the form a + A/rj in the fitting region roughly 0.5ro < 77 < 2r$. The 
parameters a and A would be identified as the string tension and the Coulomb coefficient. The 
systematic uncertainty due to the choice of fit range is small, and at most the same size as the 



Table 1: Lattice parameters. The scale a~ l is determined from the hadronic radius ro- The 

mean-field values are in the Landau gauge. The statistical uncertainty of u T is less than the 
last digit. 

(3 1G size r a' 1 [GeV] u a u T tj M f 

5.75 3.072 12 3 x 96 2.786(15) 1.100( 6) 0.7620(2) 0.9871 1.2953(4) 

5.95 3.1586 16 3 x 128 4.110(23) 1.623( 9) 0.7917(1) 0.9891 1.2494(2) 

6.10 3.2108 20 3 x 160 5.140(32) 2.030(13) 0.8059(1) 0.9901 1.2285(2) 



statistical error. Table [| shows the value of ro and the a~ l determined by setting the physical 
value of ro as r$ = 395 MeV (ro — 0.5 fm). The quoted error represent only the statistical 
uncertainty. 

The mean-field values, u a and u T , are obtained as the average of the link variables in the 
Landau gauge, where the mean-field values are used self-consistently in the fixing condition 
||. The results are also listed in Table [l]. The mean-field value of the temporal gauge field 
has small error, and close to unity, tjmf = Ur/ua, the mean-field estimate of r\ = £/7g> 
is close to the value of r\ determined nonperturbatively by Klassen. This suggest that the 
tadpole-improvement works well also on the anisotropic lattice. 

4.2 Quark field calibration 

As described in Section pi we use the relativistic dispersion relation in the calibration of param- 
eters in the quark action. Since the gauge field calibration is in the accuracy of one percent, 
we aim to tune the quark parameters to a similar level. 

For convenience, we choose k and jf as the input parameters and determine K a and k t 
along eq. (|2.5|). Fixing k corresponds to fixing the bare quark mass in the spatial lattice unit. 
For each values of (k, 7f), the pseudoscalar and vector meson correlators are obtained with zero 
and finite momenta. The fermionic anisotropy £f is defined through the relativistic dispersion 



relation, eq. (3.2). We assume the linear dependence of £p on 7^ in the vicinity of £p ~ £. We 



use linear interpolation to obtain 7^, the value of jp with which the relation £p = £ holds. 

Result of the dispersion relation. The parameters (k,jf) used in the calibration are 
listed in Table ^, [| and ||] for j3 = 5.75, 5.95 and 6.10, respectively. As the meson operators 
at the source, we adopt the smeared operators with appropriate smearing functions. For the 
light quark region, the Gaussian function is used as the smearing function, with the width of 
0.2 ~ 0.4 fm. In the charm quark mass region, we also use the measured wave function for the 
smearing function. We measure the two point functions for momenta p = n ■ (2ir/L), where L 
is spatial lattice extent and n = (0, 0, 0), (1, 0, 0), (1, 1, 0), (1, 1, 1) and (2, 0, 0). All rotationally 
equivalent n's are averaged. The standard procedure is used in extracting the energy at each 
momentum. 

The energies are then fitted to the linear or quadratic forms in p 2 to extract the fermionic 
anisotropy £p in each channel. In the case of linear fit, we use only three lowest momentum 
states, n = (0,0,0), (1,0,0) and (1,1,0). We assume the linear dependence of £f(jf) in Jf, 
and this is indeed verified in several examples. 

Figure shows typical effective mass plots for pseudoscalar and the vector mesons. The 
energies of finite momentum states are successfully extracted from the region in which the 
correlator shows plateaus except the lightest quark region, k < 0.120, which severely suffer 




Figure 2: Effective mass plots of PS and V mesons at K = 0.105, jp = 3.8 on (3 = 5.75 
lattice. From bottom to top, the states with integer momentum vectors n = (0,0,0), (1,0,0), 
(1, 1, 0), (1, 1, 1) and (2, 0, 0). Horizontal solid lines represent the result of fit of correlators and 
the fitting range. The statistical errors are represented by the dashed lines. The state with 
n = (2, 0, 0) state is not used in the fit of the dispersion relation at this /3. 
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Figure 3: Dispersion relations of PS and V mesons. The left figure shows the data at k = 0.105, 
7_F = 3.8 and 3.7 on (5 = 5.75 lattice. The right figure shows those at k = 0.115, 7_f = 4.0 
and 4.1 on (5 = 6.10 lattice. Solid lines represent the linear fit and the dashed lines show the 
quadratic fit. 



from the statistical fluctuation. 

The dispersion relation for k = 0.105 at /3 = 5.75 is shown in the left figure of Figure [3|. 
Because of rather large lattice artifact, the fit to the quadratic form in p 2 with the energy 
of n = (2,0,0) state is not a good description of data. We therefore use only four lowest 
energy states in the quadratic fit. On the other hand, the results of the linear fit (with lowest 
three states) and the quadratic fit coincides with good accuracy. For a few largest hopping 
parameters, higher momentum states suffers so large statistical fluctuations that we inevitably 
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Figure 4: 7^ dependence of £p (left) and meson masses (right) at k = 0.105 on (3 = 5.75 
lattice. 

adopt the linear fit. Since at other values of k the resultant £p coincides with quadratic fit, 
we adopt the linear fit for all values of k at this (3. Since the same method as (3 = 5.75 case is 
adopted at (3 = 5.95, we do not repeat the explanation for the fitting procedure. 

The right figure of Figure |3| shows the dispersion relation of mesons at (3 = 6.10 and 
k = 0.115, which corresponds to a similar quark mass as k = 0.105 at (3 = 5.75. The dispersion 
relation is much improved, and the quadratic fit is successfully applied including n = (2, 0, 0) 
state. Although the difference between the linear fit is small, as is shown in the figure, we adopt 
the result of quadratic fit to determine £p except for most light quark region. For these three 
largest k, correlators with n = (2, 0, 0), and occasionally (1, 1, 1), suffer from so large statistical 
fluctuations that the energy of the states cannot be reliably extracted. In these cases, we fit 
the energy to the linear form. 



Calibration of 7^ for each quark mass. In the left figure of Figure ||, £p is plotted as 
the function of ■jf for k = 0.105 at (3 = 5.75. It is clear that £p depends linearly on jp. The 
results at k = 0.101, 0.110 and 0.112 also show similar behavior. We therefore assume the 
linear dependence also for other values of k, and interpolate £p to find out 7^ in each channel. 
The numerical results of 7^ , 7^ and 7^ as the average over PS and V mesons are listed 

in Table 0-f§ for each (3. These tables also show the interpolated masses of PS and V mesons. 

*(PS) *(V) 

We find a tendency that 7^ is slightly larger than 7^ in whole k region. This deviation 

seems to become smaller for larger (3. The reason for the discrepancy is understood as the 

systematic errors of 0(aa), which will be examined in the next subsection in detail. 

We also plot the 7^ dependence of meson masses at k = 0.105 and (3 = 5.75 in the right 

figure of Figure ||[ This shows that the meson mass is linear in 7^ in this range, and linear 

interpolation is successfully applicable to determine the meson masses at 7^. Although the 

jp dependence of meson mass is in general unknown for other region of k, we expect that 

the linear interpolation would work in good accuracy. The meson masses at 7^ are also listed 

in Table |2|-|| for each (3. How the uncertainty in 7^ affect on the spectrum is an important 

problem, which will be examined in the next subsection. 



Table 2: Calibration parameters and result at j3 = 5.75. Linear fit is applied to the dispersion 
relation in the determination of £p. 



input jp 



N, 



conf 'Yp 



WW 



If 



WT 



r F 



mpsilp) m v(l F ) 



0.124 
0.122 
0.120 
0.118 
0.116 
0.114 
0.112 
0.110 
0.105 
0.101 
0.097 
0.095 
0.093 



3.9, 4.0 
3.9, 4.0 
3.9, 4.0 
3.9, 4.0 
3.9, 4.0 
3.9, 4.0 
3.8,3.9,4.0 
3.8,3.9,4.0 

3.7,3.8,3.9,4. 

3.7,3.8,3.9,4. 
3.5, 3.6 
3.5, 3.6 
3.5, 3.6 



400 
400 
400 
400 
300 
200 
200 
200 
160 
160 
160 
160 
160 



3.935(77) 
3.904(48) 
3.892(43) 
3.906(36) 
3.875(35) 
3.899(36) 
3.854(29) 
3.878(33) 
3.807(30) 
3.730(26) 
3.626(19) 
3.579(18) 
3.530(17) 



3.83(18) 
3.884(82) 
3.888(54) 
3.894(42) 
3.841(42) 
3.842(47) 
3.806(37) 
3.827(41) 
3.754(37) 
3.679(31) 
3.587(24) 
3.540(23) 
3.490(21) 



3.919(72) 
3.899(45) 
3.891(38) 
3.901(31) 
3.861(33) 
3.878(36) 
3.836(30) 
3.857(34) 
3.786(31) 
3.709(27) 
3.611(20) 
3.564(19) 
3.514(18) 



0.1497( 6) 
0.2044( 4) 
0.2523( 8) 
0.2967( 9) 
0.3408(13) 
0.3819(17) 
0.4252(17) 
0.4654(22) 
0.5738(28) 
0.6653(30) 
0.7647(28) 
0.8166(29) 
0.8704(29) 



0.2294(17) 
0.2650(12) 
0.3018(12) 
0.3387(12) 
0.3774(15) 
0.4142(19) 
0.4546(18) 
0.4918(23) 
0.5954(28) 
0.6845(30) 
0.7823(28) 
0.8333(29) 
0.8864(29) 



Table 3: Calibration parameters and result at (3 = 5.95. Linear fit is applied to the dispersion 
relation in the determination of £p. 



K 


input jf 


N t 

1 v con j 


*(PS) 

7f 


*(V) 

7f 


If 


mps{l F ) 


mvilp) 


0.124 


3.9, 4.0 


500 


4.073(95) 


4.15(12) 


4.103(80) 


0.1177( 6) 


0.1649( 9) 


0.123 


3.9 


4.0 


500 


4.041(69) 


4.095(93) 


4.060(60) 


0.1456( 2) 


0.1847( 8) 


0.122 


3.9 


4.0 


500 


4.029(55) 


4.076(68) 


4.048(48) 


0.1712( 4) 


0.2045( 8) 


0.120 


3.9 


4.0 


500 


4.019(36) 


3.996(55) 


4.012(34) 


0.2186( 6) 


0.2444( 9) 


0.118 


3.9 


4.0 


500 


4.003(29) 


4.000(35) 


4.002(28) 


0.2625( 8) 


0.2841( 9) 


0.115 


3.9 


4.0 


360 


3.992(28) 


3.969(36) 


3.983(29) 


0.3260(12) 


0.3431(13) 


0.110 


3.9 


4.0 


300 


3.945(28) 


3.946(36) 


3.946(29) 


0.4297(19) 


0.4427(19) 


0.107 


3.9 


4.0 


300 


3.910(25) 


3.911(31) 


3.910(26) 


0.4930(20) 


0.5046(21) 


0.104 


3.9 


4.0 


200 


3.876(28) 


3.875(36) 


3.876(30) 


0.5573(27) 


0.5677(27) 


0.102 


3.9 


4.0 


200 


3.848(26) 


3.847(34) 


3.847(28) 


0.6016(28) 


0.6113(28) 


0.100 


3.9 


4.0 


200 


3.815(25) 


3.816(32) 


3.815(27) 


0.6470(29) 


0.6562(29) 


0.097 


3.8 


3.9 


200 


3.766(24) 


3.765(30) 


3.766(26) 


0.7178(32) 


0.7264(32) 


0.093 


3.7 


3.8 


200 


3.688(23) 


3.687(29) 


3.688(25) 


0.8180(37) 


0.8257(37) 



Fit of 7£. To represent 7^ as the function of re, we introduce the quark mass as 



1 /l 



m n 



2£ \k 



(4.1) 



This is similar relation as mo, the bare quark mass in the temporal lattice unit, while the 
present form is independent of 7^. re c is determined from massless point of the pseudoscalar 
meson mass. We extrapolate mp s linearly in 1/re using two largest values of re, and find as 
re c = 0.12640(5) at /3 = 5.75, 0.12592(6) at (3 = 5.95 and 0.12558(4) at = 6.10. 

In the calibration at each re, we found that the value of 7^ is easily determined precisely (to 
the level of 1 %), while it become difficult with increasing re toward re c . However, it is expected 
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Table 4: Calibration parameters and result at j3 = 6.10. For k =0.124, 0.123 and 0.122, 
the dispersion relation is fitted to the linear form in determining £p. For the remaining n's, 
quadratic fit is applied. 



K 


input jp 


N t 

1 y con j 


*(PS) 

7f 


*{V) 

7f 


If 


m P s{l F ) 


mvil F ) 


0.124 


4.0, 4.1 


600 


4.020(77) 


3.63(27) 


3.991(79) 


0.1008( 5) 


0.1379( 7) 


0.123 


4.0 


4.1 


600 


4.005(52) 


3.860(98) 


3.973(55) 


0.1294( 1) 


0.1582( 5) 


0.122 


4.0 


4.1 


600 


3.998(41) 


3.913(69) 


3.976(44) 


0.1549( 3) 


0.1787( 6) 


0.120 


4.0 


4.1 


400 


4.040(29) 


4.064(45) 


4.047(30) 


0.2007( 6) 


0.2182( 7) 


0.118 


4.0 


4.1 


200 


4.040(33) 


4.014(52) 


4.032(34) 


0.2440( 9) 


0.2579(10) 


0.115 


4.0 


4.1 


200 


4.024(28) 


4.008(41) 


4.019(30) 


0.3067(12) 


0.3176(13) 


0.110 


4.0 


4.1 


200 


4.013(38) 


3.996(54) 


4.007(43) 


0.4078(26) 


0.4160(27) 


0.107 


4.0 


4.1 


200 


3.988(36) 


3.986(46) 


3.987(39) 


0.4694(29) 


0.4766(30) 


0.104 


4.0 


4.1 


200 


3.951(33) 


3.955(42) 


3.952(36) 


0.5331(31) 


0.5396(31) 


0.102 


3.9 


4.0 


200 


3.918(32) 


3.928(40) 


3.922(35) 


0.5773(34) 


0.5834(34) 


0.100 


3.9 


4.0 


200 


3.877(25) 


3.883(32) 


3.879(27) 


0.6240(28) 


0.6298(28) 


0.097 


3.9 


4.0 


200 


3.834(22) 


3.839(28) 


3.836(24) 


0.6931(28) 


0.6984(28) 


0.093 


3.8 


3.9 


200 


3.769(21) 


3.776(26) 


3.772(23) 


0.7903(32) 


0.7953(31) 



7^ smoothly approaches to certain definite value in the limit of m q — ► 0, since in this limit our 
form of action is simply a direct generalization of clover quark action on the anisotropic lattice. 
In fact, as shown in Subsection [T^, 7^ linearly depends on itlq at the tree level. In taking 
the limit of m q — » for j F , the precise values of mean-filed values do not matter, since in 
the definition of -yp, tadpole-improvement is applied as a multiplicative factor, t\mf = u T /u a . 
Therefore the most reliable way to determine the value of 7^ in the light quark region is a 
global fit of j F , assuming appropriate form of m q dependence. 

We fit the result of I/7J? to the linear form in m q and the quadratic form in m q . The result 
of fit is listed in Table || and also shown in Figures || as the solid and the dashed curves. Since 
the obtained points for 7^ is with different number of configurations, these points correlate 
in not obvious way. We quote the errors and x 2 °f uncorrelated fit in Table ISL As shown 
in these table and figures, the linear fit in rrii well represents the data. The values of £o> 
which is the 1/yp in the chiral limit, is close to the tree level result l/£ = 1/4, which implies 
that the tadpole-improvement works well. Apart from the lightest quark mass region, 7^ is 
determined within 1 % accuracy, and there the curves of two fits are consistent with each other. 
In approaching the chiral limit, there is a systematic error concerning the fit form, as well as 
the statistical error. We estimate the latter by the error of fit in £0 (from the quadratic fit in 
m q ), and as about 1 %. This relatively small statistical error is due to the global fit of 1/yp 
with assumed form of m q dependence. The systematic error in adopting specific form of fit is 
estimated by the difference between these two fits, and also results in 1 % level. Adopting the 
linear form in m q , we conclude that the 7^ is determined under the assumed dispersion relation 
within 1 % statistical accuracy in the whole quark mass region less than around the charm 
quark mass, while in the chiral limit there is additional 1 % systematic uncertainty concerning 
the form of fit. 
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P=5.75 
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P=5.95 

o data 
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--- fit: y f * <o + ^ m q <2 m 
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p=6.10 
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— fit:y F "^o + ^ m q 2 

-- fit: y F *~ =£„+£, m q +C 2 m q ' 
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1/k 



Figure 5: l/^fp vs 1/k at each (3. Solid lines show the fit linear in mi while the dashed lines 
represent the fit quadratic in m q . 
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Table 5: Fit results of 7^. 






fit type 


Co 


Ci 


C2 


X 2 /N df 


l* F {m q = 0) 


5.75 


linear 
quad. 


0.2558( 9) 
0.2564(23) 


-0.007(28) 


0.230(12) 
0.247(68) 


1.83 / 11 
1.77 / 10 


3.909(14) 
3.901(34) 


5.95 


linear 
quad. 


0.2490( 8) 
0.2465(18) 


0.036(23) 


0.189(15) 
0.095(61) 


3.52 / 11 
1.01 / 10 


4.016(13) 
4.057(30) 


6.10 


linear 
quad. 


0.2479( 9) 
0.2493(18) 


-0.022(24) 


0.143(14) 
0.200(63) 


4.44 / 11 
3.55 / 10 


4.034(14) 
4.011(28) 



4.3 Uncertainties in calibration 

In the last subsection, we determined 7^ as a global function of m q . This expression inevitably 
suffer from systematic uncertainties as well as the statistical uncertainty: 



If 



7 *(prop) + ^(stat) + ^(O(aa)) + ^ 



(,0(a 2 )) f ,,(chiral), 



'(+H 



(4.2) 



7^ represent the proper value of the bare anisotropy. 5^/ F is the statistical error in 

determination of 7^, and at 1 % level. The last two terms are main sources of systematic errors 
due to finite lattice artifact. The first one, Sjp , is from the tree-level approximation of 

clover coefficients. We estimate the size of this error by the difference between the values of 
j F determined with PS and V mesons. The second systematic error, 5^/ F , is estimated 
by comparing the results of calibration from two different forms of dispersion relations which 
differ by 0(a 2 ). In addition to these systematic uncertainties, in the chiral limit there is also 
a systematic error concerning the form of fit of 7^ in m q . 

Another important subject is to estimate how the observables are affected by the uncertainty 
in 7^. We study the response of meson masses with respect to the change of 7^ at each k, 
from which the effect of 7^ on meson masses for given quark mass is approximately estimated. 
Strictly speaking, changing 7^ for a fixed k induces a slight change in the quark mass, hence 
the above analysis is adequate for relatively heavier quark mass region, such as m s < m q . We 
postpone the study of the effect on the light hadron spectrum to the end of next section. 



Difference between 7^ ' and 7^ . Since we use the 0(a) improved quark action, the 
main contribution from the 0(a) lattice artifact is absent. However, since the clover coef- 
ficients are not tuned beyond the tree-level, 0(aa) error still remains, although the tadpole- 
improvement would partially removes this effect. An appropriate probe of this systematic effect 

on the calibration is the difference between 7^'s of the pseudoscalar and the vector mesons. 

n *(PS) *(V} 

Figure |6| shows 5j F = 7^ — 7 F . At p = 5.75, there is a systematic difference of 5"f F 

from zero except the small quark mass region, where the statistical error is dominant. At 

j3 = 5.95, 5j F is consistent with zero in whole k region. This implies that the O(aa) error in 

the calibration is sufficiently reduced at this (3. In the case of (3 = 6.10, ^7^ is also consistent 

with zero except the lightest quark region. In this region, precise determination of energy at 

finite momentum states is difficult due to the statistical fluctuation and hence the resultant jf 

contains large uncertainty. We regard that the O(aa) effect in the calibration is also small at 

this (3. 
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1/K 

*(PS) *(V) 

Figure 6: 5~f F = ^ F — ^ F at =5.75 and 5.95 in upper part, and 6.10 in lower part. 

0(a?) systematic uncertainty. Although we employed the continuum dispersion relation, 
this introduces a systematic error of 0(a 2 ) to the calibration of jp. In order to estimate the 
typical size of this error, we compare 7^ determined above with 7p( KG \, the result obtained 
using the dispersion relation from the lattice Klein-Gordon action, ( p.3[ ). Figure ffi shows this 
comparison at (3 = 5.75 and 6.10. In extracting £,kg from the Klein-Gordon dispersion relation, 
we fit cosYiE(p) to the linear form in X^sm 2 (pi/2) using lowest three momentum states. As 
explained in Section ||, expected difference of £p and £kg is 0(m 2 ), where m is the meson 
mass. Although the explicit relation between ^pi^o and 7^ is unknown, one can expect that 
the difference of them is also in 0(m 2 ), and hence increase with increasing quark mass. This 
behavior is clearly observed in Figure 0. Table ^ is the result of fit of Ipt^a) ^° ^ rie nnear fo rm 
in m 2 

We find small difference between the results with relativistic and Klein-Gordon dispersion 
relations in the small quark mass region. This difference decreases as increasing /?, and seems 
to be sufficiently reduced at (3 = 6.10. The typical size of difference in 7p's at the chiral limit 
is less than 3, 2, 1 % at j3 = 5.75, 5.95 and 6.10, respectively. The important feature is that 
two procedures tend to coincide with each other with increasing (5. We also observe that the 
difference between 7^ and ^pi KG \ increases in the large quark mass region, m q > 0.2a" 1 . This 
is consistent behavior, since there the Klein-Gordon dispersion relation fails to incorporate 
the quark mass dependence properly, and ^f*pi KG \ is expected to be larger than 7^, in 0(m 2 ). 
Therefore we conclude that the uncertainty due to the assumed form of the meson dispersion 
relation is under control and smoothly disappears in approaching the continuum limit. 

Uncertainty in meson mass due to calibration error. Another important issue on 
the systematic errors is how the uncertainty in 7^ is transmitted to the observables. As an 
important example, here we focus on the effect on the meson masses. 

Since we linearly interpolate the meson masses in 7^, we obtain dm/d'jF at 7^ = 7^ from 
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Table 6: The result of linear fit in nil of Jp/^Q), the tuned bare anisotropy with the Klein- 
Gordon dispersion relation. 






fit type 


Co 


c 2 


X Z /N df 


l* F {m q = 0) 


5.75 


linear 


0.2488( 8) 


0.112(11) 


2.33 / 11 


4.019(13) 


5.95 


linear 


0.2446( 8) 


0.071(14) 


2.14 / 11 


4.088(13) 


6.10 


linear 


0.2484(10) 


0.039(15) 


1.59 / 11 


4.026(16) 





P=5.75 






, 


0.28 


O Relativistic DR 

■ Lattice Klein-Gordon 






- 


0.27 








- 


0.26 
0.25 


- |&¥r r 


-r 


A* 


i'' 



0.28 



0.27 



0.26 



0.25 



0.24 



0.23 



[3-6.10 

O Relativistic DR 

■ Lattice Klein-Gordon 




■■■■■■ 



7.5 8.0 8.5 9.0 9.5 10.0 10.5 11.0 11.5 7.5 8.0 8.5 9.0 9.5 10.0 10.5 11.0 11.5 

1/k 1/k 

Figure 7: Comparison of the results of calibrations with two types of dispersion relations. The 
curves represent the the result of linear fit in ?tt.q. Solid lines are the fit results with relativistic 
dispersion relation and the dashed lines are those with Lattice Klein-Gordon dispersion relation. 

the slope of the linear fit. In Figure |8], we show dm/d'jp, the response of meson mass to the 
bare anisotropy, at (3 = 6.10 in two ways. Similar feature is found in the results at = 5.75 
and (3 = 5.95. In the left figure dm/djp is shown as the function of 1/k. In the case of vector 
meson, it seems to decrease linearly in quark mass from zero at the massless limit. On the 
other hand, for the pseudoscalar meson, dm/djp slightly positive in the vicinity of l/n c . This 
behavior may be due to the uncertainty in the definition of k, because if k is properly related to 
the fixed quark mass, increasing jp implies to increase the propagation in temporal direction, 
hence it corresponds to decreasing quark mass. Therefore, the present analysis may not be 
adequate for estimating the response of masses with respect to jp in the vicinity of the chiral 
limit. Observing Figure ||, one can find that the range of quark mass larger than the strange 
quark mass do not suffer from the ambiguity in the definition of n. 

We have no clear explanation for that dmjdjp seems to be proportional to the quark mass. 
Practically, it is a good feature that the ambiguity oijp has only little effects on the meson mass 
in the small but nonzero quark mass region, since there relative change of mass is significant. 
Except the lightest quark mass region, the determination of j F is directly performed with 1 
% accuracy, which means the uncertainty of 7^ is around 0.04. The right figure in Figure |8| 
implies that the uncertainty in the meson mass is less than 1 %. This uncertainty is most 
severe in the heavy quark region, and become milder as quark mass decreases. 

While we find that the meson masses at certain k is not so sensitive to the uncertainty 
of jp, the same argument does not hold for the chiral limit. Since the pseudoscalar meson 
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Figure 8: The response of meson masses to the change of jp at (3 = 6.10. The left figure 
shows dm/d'jF, while the ratio (l/m)dm/djF is shown in the right figure, the result for the 
pseudoscalar and the vector mesons are represented by circles and triangles respectively. 

mass becomes zero in the chiral limit, the relative uncertainties 5mps/mps for a fixed k near 
the chiral limit is of course very large. However, this is not the correct way to estimate the 
uncertainties in the mass spectrum in the chiral limit. What one is really interested in for the 
chiral limit is not the change the hadron masses including the pion mass for a fixed n but the 
change of hadron masses except the pion mass at the point where the pion becomes massless. 
Since the critical hopping parameter k c is affected by the change in jp, one need to treat the 
chiral limit carefully. In the next section, we discuss the uncertainties of the hadron spectrum 
in the chiral limit based on the extrapolation in terms of the pseudoscalar meson mass squared 
instead of 1/k. 



4.4 Summary of calibration 

In this section, we have implemented the anisotropic 0(a) improved Wilson action in the 
region of quark mass up to around the charm quark mass, at three values of (3 at £ = 4.0. 
The fermionic anisotropy £p is extracted from the meson dispersion relation. Then we find the 
value of bare anisotropy parameter, y F , at which £p = £ holds. The value of 7^ in the massless 
limit is obtained by extrapolating the data by fitting to the linear form in mz, where m q is 
naively defined quark mass. This is the most reliable way to determine 7^ for small quark 
mass region, since there the statistical fluctuation in finite momentum states is severely large. 
The fit of l/7jp to the linear form in nil seems quite successful, and 7^ at the chiral limit is 
close to the tree-level value, £. The statistical uncertainty in 7^ is estimated as in the order 
of 1 % in whole explored quark mass region. In the chiral limit, there is also 1 % systematic 
uncertainty concerning the form of fit. Here we summarize the main result of calibration, the 
expression of 7^ for a given k: 

1 , , o 1/11 



If 2 ? 

P = 5.75 : Co = 0.2558(9), (2 = 0.230(12), 

(3 = 5.95 : Co = 0.2490(8), C2 = 0.189(15), 

P = 6.10 : Co = 0.2479(9), C2 = 0.143(14), 



0.12640(5) 
0.12592(6) 
0.12558(4). 



(4.3) 

(4.4) 
(4.5) 
(4.6) 
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Table 7: Quark parameters used in the hadron spectroscopy. 

j5 7_F values of k -^conf 

5.75 3.909 0.1240, 0.1230, 0.1220, 0.1210 200 

5.95 4.016 0.1245, 0.1240, 0.1235, 0.1230 100 

6.10 4.034 0.1245, 0.1240, 0.1235, 0.1230 100 



To examine the uncertainty in the calibration, we have also carried out the following anal- 
ysis, (i). The difference between 7^ for the pseudoscalar and vector mesons, which signals the 
0(aa) systematic error. We observed that this difference decreases as decreasing lattice spac- 
ing, and already consistent with zero at (3 = 5.95. (ii). Comparison of 7^ with the continuum 
and the Klein-Gordon dispersion relations. This is for a estimate of the size of 0(a 2 ) system- 
atic uncertainty. The results with two dispersion relations tends to coincide with each other 
in decreasing the lattice spacing. The behavior in the large quark mass region is consistent 
with the expected behavior, (iii). Response of meson masses to the change of jf- The effect 
of uncertainty of 7^ on the meson masses is less than 1 %, if 7^ is determine at this accuracy. 
This result is applicable to the relatively heavier quark mass region, such as m s < m„, and 
therefore in the region, the errors in the calibration is under control. 

5 Light hadron spectroscopy 

In this section, we apply the result of the last section to the calculation of light hadron spectrum. 
Our analysis is performed in two steps: 

(i) By taking the central value of jp = 7^ we obtain the light hadron masses in the strange 
quark mass region m s < m q < 2m s . By extrapolating masses in mpg, the hadron spectrum 
at the physical light quark masses are determined. We compare our result with the result 
by UKQCD Collaboration fllql , which has been obtained on the isotropic lattice with O(a) 
improved quark action. 

(ii) We study the response of the light hadron spectrum to the change of the anisotropic pa- 
rameter 7jp — » 7jp + 57_F- The extrapolation in mp s is significant to circumvent the uncertainty 
in the definition of n. 

5.1 Calculation of hadron spectrum 

The spectroscopy of light hadrons are performed on the same lattices used in the calibration, 
while with smaller numbers of configurations. The parameters are listed in Table [7|. At each 
/3, we use four values of k corresponding to the quark mass of m s ~2m s . In this region, we 
regard that m q is sufficiently small so that we can adopt the value of "f F in the massless limit. 
Therefore the bare anisotropy is set to the central value of 7^ at m q = 0, which is determined 
in the calibration as the linear form in m?. 

We use the standard hadron operators and procedure to extract the hadron masses. The 
quark propagators are smeared at the source with Gaussian smearing function with the devia- 
tion ~ 0.4 fm, in the Coulomb gauge. The periodic boundary condition is adopted in all four 
directions for the quark field. For baryons, two of quarks are treated as they has degenerate 
masses. Then the quark content of baryons is specified by two k's, k\ and K2, for a pair of 
quarks and the other quark, respectively. Figure |9| shows the effective mass plot for octet and 
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Figure 9: Effective mass plots for octet and decuplet baryon correlators with degenerate quark 
masses at /3 = 5.75. Horizontal solid line represent the fit range and the mass from the single 
exponential fit. 



Table 8: Hadron spectrum at (3 = 5.75. When quark masses are degenerate, i.e. K\ 
X-type and the A-type octet baryon correlators are identical. 



k,2, the 



«1 


K2 


m PS 


m v 


m oei(£) 


m oct(A) 


m d ec 


0.1210 


0.1210 


0.22909(47) 


0.2821(10) 


0.4257(17) 


- 


0.4564(24) 


0.1210 


0.1220 


0.21716(48) 


0.2730(11) 


0.4146(18) 


0.4161(18) 


0.4472(25) 


0.1210 


0.1230 


0.20501(51) 


0.2640(12) 


0.4034(19) 


0.4068(19) 


0.4383(27) 


0.1210 


0.1240 


0.19260(54) 


0.2552(13) 


0.3920(20) 


0.3977(21) 


0.4299(29) 


0.1220 


0.1210 


- 


- 


0.4059(19) 


0.4042(19) 


0.4381(27) 


0.1220 


0.1220 


0.20480(50) 


0.2637(12) 


0.3946(19) 


- 


0.4289(28) 


0.1220 


0.1230 


0.19214(52) 


0.2546(13) 


0.3831(20) 


0.3852(21) 


0.4199(30) 


0.1220 


0.1240 


0.17911(55) 


0.2458(15) 


0.3713(21) 


0.3760(22) 


0.4114(33) 


0.1230 


0.1210 


- 


- 


0.3863(21) 


0.3820(20) 


0.4202(30) 


0.1230 


0.1220 


- 


- 


0.3747(22) 


0.3724(21) 


0.4109(32) 


0.1230 


0.1230 


0.17886(54) 


0.2454(15) 


0.3629(23) 


- 


0.4018(35) 


0.1230 


0.1240 


0.16503(57) 


0.2364(17) 


0.3506(24) 


0.3538(25) 


0.3932(39) 


0.1240 


0.1210 


- 


- 


0.3674(25) 


0.3587(22) 


0.4032(37) 


0.1240 


0.1220 


- 


- 


0.3555(26) 


0.3490(24) 


0.3938(40) 


0.1240 


0.1230 


- 


- 


0.3432(27) 


0.3394(25) 


0.3845(44) 


0.1240 


0.1240 


0.15015(60) 


0.2272(21) 


0.3303(28) 


- 


0.3757(51) 



decuplet baryons with degenerate quark mass case, k\ = Ki- The meson correlators are fitted 
to the single hyperbolic cosine form. For baryons, we apply the single exponential fit in the 
region in which there is negligible contribution of the negative parity baron from the other 
temporal boundary. The result of fit is listed in Tables ||- [[(]. For mesons, since the order of 
K\ and K2 is unimportant, the masses at the exchanged set of (k\, K2) are omitted. The masses 
of A-type octet baryon at degenerate (k\, K2) are also omitted, since they are identical to the 
masses of S-type. 
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Table 9: Hadron spectrum at (3 = 5.85. 



«1 


K2 


m PS 


my 


m oct(T,) 


m oct(A) 


m d ec 


0.1230 


0.1230 


0.14580(40) 


0.1853(10) 


0.2788(17) 


- 


0.3036(25) 


0.1230 


0.1235 


0.13896(41) 


0.1804(11) 


0.2726(17) 


0.2734(18) 


0.2987(26) 


0.1230 


0.1240 


0.13196(43) 


0.1755(11) 


0.2662(18) 


0.2680(19) 


0.2938(27) 


0.1230 


0.1245 


0.12477(46) 


0.1707(13) 


0.2596(19) 


0.2625(20) 


0.2890(30) 


0.1235 


0.1230 


- 


- 


0.2676(18) 


0.2666(18) 


0.2937(27) 


0.1235 


0.1235 


0.13190(43) 


0.1754(11) 


0.2612(19) 


- 


0.2888(29) 


0.1235 


0.1240 


0.12462(45) 


0.1705(12) 


0.2546(19) 


0.2557(20) 


0.2839(31) 


0.1235 


0.1245 


0.11709(47) 


0.1656(14) 


0.2478(20) 


0.2501(21) 


0.2790(33) 


0.1240 


0.1230 


- 


- 


0.2562(20) 


0.2539(19) 


0.2839(31) 


0.1240 


0.1235 


- 


- 


0.2496(21) 


0.2484(20) 


0.2789(33) 


0.1240 


0.1240 


0.11702(47) 


0.1654(14) 


0.2427(21) 


- 


0.2740(35) 


0.1240 


0.1245 


0.10908(50) 


0.1605(16) 


0.2356(22) 


0.2370(23) 


0.2692(38) 


0.1245 


0.1230 


- 


- 


0.2445(23) 


0.2405(21) 


0.2743(37) 


0.1245 


0.1235 


- 


- 


0.2376(24) 


0.2348(22) 


0.2693(39) 


0.1245 


0.1240 


- 


- 


0.2305(25) 


0.2289(24) 


0.2644(42) 


0.1245 


0.1245 


0.10063(53) 


0.1555(19) 


0.2229(26) 


- 


0.2597(47) 



Table 10: Hadron spectrum at (3 = 6.10. 



«1 


K2 


mps 


my 


"W(E) 


™>oct(k) 


m d ec 


0.1230 


0.1230 


0.12950(29) 


0.1587( 6) 


0.2394(11) 


- 


0.2590(17) 


0.1230 


0.1235 


0.12284(30) 


0.1538( 6) 


0.2332(12) 


0.2340(12) 


0.2542(18) 


0.1230 


0.1240 


0.11603(31) 


0.1491( 7) 


0.2269(12) 


0.2288(12) 


0.2495(19) 


0.1230 


0.1245 


0.10904(33) 


0.1446( 8) 


0.2204(13) 


0.2236(13) 


0.2451(21) 


0.1235 


0.1230 


- 


- 


0.2283(12) 


0.2273(12) 


0.2493(19) 


0.1235 


0.1235 


0.11595(31) 


0.1489( 7) 


0.2219(13) 


- 


0.2445(20) 


0.1235 


0.1240 


0.10886(32) 


0.1442( 8) 


0.2154(13) 


0.2166(13) 


0.2399(22) 


0.1235 


0.1245 


0.10153(33) 


0.1398( 9) 


0.2087(14) 


0.2112(14) 


0.2355(24) 


0.1240 


0.1230 


- 


- 


0.2171(14) 


0.2148(13) 


0.2400(22) 


0.1240 


0.1235 


- 


- 


0.2106(14) 


0.2093(14) 


0.2352(24) 


0.1240 


0.1240 


0.10142(33) 


0.1395( 8) 


0.2038(14) 


- 


0.2306(26) 


0.1240 


0.1245 


0.09366(35) 


0.1352(10) 


0.1968(15) 


0.1983(16) 


0.2263(30) 


0.1245 


0.1230 


- 


- 


0.2058(16) 


0.2017(15) 


0.2314(28) 


0.1245 


0.1235 


- 


- 


0.1990(16) 


0.1960(15) 


0.2267(30) 


0.1245 


0.1240 


- 


- 


0.1919(17) 


0.1903(16) 


0.2221(34) 


0.1245 


0.1245 


0.08535(36) 


0.1310(12) 


0.1845(18) 


- 


0.2180(41) 



5.2 Extrapolation to the chiral limit 

In order to avoid the ambiguities in the definition of the quark mass, we extrapolate the hadron 
masses to the chiral limit in terms of the pseudoscalar meson mass squared, instead of 1/re. 
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We assume the relation 

mp S (mi,m 2 )=B'(mi+ra 2 ), (5.1) 

then for the degenerate quark masses, mi = 771,2, m ps = 2-Brai holds. Then instead of 
777j (i=l,2), one can use mps(mi,mi) 2 as the variable in the chiral extrapolation. For other 
hadrons, vector meson and octet and decuplet baryons, we also use the linear relations: 

777^(7771,7712) = mv(0,0) + By ■ (mi +m 2 ), (5.2) 

m oct (mi,m2,m 3 ) = m oct (0,0,0) + B oct ■ (mi + m 2 + m 3 ), (5.3) 

■m de c('mi,m2,m 3 ) = m dec (0, 0,0) + B dec ■ (mi + m 2 + m 3 ). (5.4) 

The hadron spectrum and the result of fit are shown in Figure |l^. The vertical axis is the 
averaged pseudoscalar meson mass squared, 



N q 



,2 



{m- ps (mi)) = — Y, m PS (mi,mi) = — 2J 25m; (5.5) 

1 i=l 5 i=l 

with N q = 2 for mesons and N q = 3 for baryons. The A-type baryon is not shown in the figure 
to avoid that the figure become too messy. The linear fit seems to be successful. 

5.3 Spectrum at physical quark masses 

To determine the hadron masses at the physical u, d and s quark mass, one need to set the scale 
of lattice. We do not distinguish the u and d quark masses, and express their mass as m n . We 
adopt two definitions, through the hadronic radius 7*0, and through the K* meson mass. These 
two methods are also adopted by UKQCD Collaboration in fi"5| , and convenient for comparison 
of our data with theirs. In |lq| , they employed the values of clover coefficient determined in 
two ways: with nonperturbative renormalization technique (NP) fl6|| , and tadpole-improvement 
(TAD) ]l2|. Then they extrapolate the masses to the continuum limit by simultaneous fit of 
these two types of data to the linear form in a 2 for NP and quadratic form in a for TAD. We 
compare our hadron spectrum at the physical quark masses with the result in the continuum 



limit of 1 15], although we ourselves do not perform the continuum extrapolation because of 



lacks of sufficient number of j3 as well as the statistical accuracy. 

Scale set by tq. The hadronic radius r^ has been already obtained in Section |3]. The 
corresponding values of spatial lattice cutoff are found in Table |]. The PS meson masses 
squared correspond to m n and m s are then defined with 777.^ = 139.6 MeV and 777^- = 495.7 
MeV (isospin averaged), respectively. These definitions are in accord with |]l^]. The hadron 
masses extrapolated or interpolated to the physical points are shown in Figure ITO and listed in 



Table 11. For comparison with the result in [15 1 , we also list the hadron masses multiplied by 



ro£ in Table 11. In the latter case, ^ appears to multiply the quantity in the spatial lattice unit 
(ro) to ones in the temporal lattice unit (masses). In our data, differences between the results 
at (3 = 6.10 and 5.95 rather large compared with the difference between (3 = 5.95 and 5.75. 
This would be partially due to the different a dependence of 0(aa) and 0(a 2 ) lattice artifact, 
and also due to the statistical fluctuation. Our results of the meson masses seem to approach 
towards the continuum results by UKQCD Collaboration on the isotropic lattices [|l5|]. 

Scale set by m^*. In the second case, 777^. = 893.9 MeV (isospin averaged) is used to set the 
lattice scale. First we interpolate the vector meson mass to the point that the ratio of PS and 
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Table 11: Hadron spectrum for physical quark masses with the scale set by tq. 







mass [GeV] 






m^rQ 




Ref. [|Lj] 




P = 5.75 


(3 = 5.95 


/3 = 6.10 


P = 5.75 


P = 5.95 


= 6.10 


cont. 


p 


0.832(11) 


0.846(17) 


0.895(12) 


2.105(28) 


2.141(42) 


2.265(30) 


2.35(16) 


K* 


0.9251(90) 


0.938(13) 


0.977(10) 


2.342(23) 


2.375(34) 


2.473(25) 


2.54(12) 


<t> 


1.0185(71) 


1.031(10) 


1.059( 8) 


2.578(18) 


2.610(26) 


2.680(21) 


2.729(77) 


N 


1.175(16) 


1.155(22) 


1.197(18) 


2.974(40) 


2.924(55) 


3.032(45) 


2.92(24) 


A 


1.260(13) 


1.253(19) 


1.283(16) 


3.190(34) 


3.172(47) 


3.248(40) 


3.22(20) 


£ 


1.281(14) 


1.268(19) 


1.302(16) 


3.242(35) 


3.211(49) 


3.295(41) 


3.23(19) 


7T' 


1.387(12) 


1.381(17) 


1.406(15) 


3.510(31) 


3.497(43) 


3.558(37) 


3.54(15) 


A 


1.403(28) 


1.440(41) 


1.521(41) 


3.552(72) 


3.65(10) 


3.85(10) 


3.86(37) 


T,* 


1.495(25) 


1.532(36) 


1.602(37) 


3.784(63) 


3.877(90) 


4.055(93) 


4.15(29) 


"* 


1.587(22) 


1.623(31) 


1.683(32) 


4.017(54) 


4.109(78) 


4.260(81) 


4.44(22) 


n 


1.679(18) 


1.715(26) 


1.763(28) 


4.249(46) 


4.341(66) 


4.464(70) 


4.72(17) 



Table 12: Hadron spectrum for physical quark masses with the scale set by thk*- The 
parameter J is also quoted, while it is dimensionless quantity. 



P = 5.75 



mass [GeV] 
P = 5.95 



P = 6.10 /3 = 5.75 



m/rriK* 
P = 5.95 



,9 = 6.10 



Ref. [1 

cont. 



56; 



P 

K* 



0.796(11) 
0.894(11) 
0.992(12) 



0.795(16) 
0.894(16) 
0.993(16) 



0.802(11) 
0.894(11) 
0.986(11) 



0.891(12) 
1.109(13) 



0.890(18) 
1.110(18) 



0.898(12) 0.921( 
1.103(12) l.HOd 



21- 



N 


1.125(15) 


1.087(21) 


1.075(16) 


1.259(17) 


1.216(23) 


1.202(18) 


1.14CiS) 


A 


1.217(14) 


1.194(20) 


1.175(15) 


1.361(15) 


1.335(22) 


1.315(17) 


1.29(1J) 


S 


1.236(14) 


1.207(20) 


1.191(16) 


1.382(16) 


1.351(23) 


1.332(17) 


1.29&I) 


H 


1.346(14) 


1.328(21) 


1.307(15) 


1.506(16) 


1.485(23) 


1.462(17) 


l-45(i^) 


A 


1.343(27) 


1.354(39) 


1.363(37) 


1.502(30) 


1.515(43) 


1.525(41) 


1.50C#) 


S* 


1.439(25) 


1.452(35) 


1.454(33) 


1.610(28) 


1.624(39) 


1.626(37) 


1.64(ii|) 


"* 


1.535(23) 


1.549(32) 


1.544(29) 


1.717(25) 


1.733(36) 


1.727(32) 


l-79(±ig) 


O 


1.631(21) 


1.647(30) 


1.634(25) 


1.825(24) 


1.842(33) 


1.828(28) 


1.93(1 8 7 ) 



J 0.3859(47) 0.3896(95) 0.3621(47) 



V meson masses is equal to the physical value of K* and K mesons. Then this vector meson 
mass defines the lattice scale. This results in the spatial lattice cutoffs 1.053(13), 1.525(27) and 
1.817(22) GeV at P = 5.75, 5.95 and 6.10, respectively. Then the values of rrip^s correspond 
to the (u,d) and s quark masses are determined with experimental K and n meson masses. The 



hadron masses at physical quark masses are listed in Table 12. We observe similar tendency 
to the case of the scale set by r$. No signal of inconsistency with the result on isotropic lattice 
is found. 
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Figure 10: The masses of vector meson and the octet and decuplet baryons together with the 
result of linear fits. Only the X-type octet baryon is shown. The filled symbols corresponding 
to the masses at the physical u, d and s quark masses with the scale set by rn. 



22 



0.5 - 



E 

"D 
> 

E 
ii 



: > 0.4 



0.3 



I 



0.0 



II 

I 



^ experiment 

O this work 

■ UKQCD(NP) 

□ UKQCD(TAD) 



f 



0.5 
a„ [GeV 



1.0 



Figure 11: The parameter J. The a a is set by using rn. Data of UKQCD Collaboration 



(square symbols) are taken from Ref. [15] on the isotropic lattices, and slightly horizontally 
shifted. 



5.4 J-parameter 

The parameter J was introduce to probe the quenching effect in 



and defined as 



J = my 



dm\ 



dm 2 PS 



(5.6) 



mv/mps=m K * /ma 



It is known that the quenched lattice simulation does not reproduce the experimental value, 



J = 0.48(2), and gives about 20 % smaller value. We show our result of J in Figure 11, as the 
function of lattice spacing determined by rrj. We find that our results are consistent to those 
by UKQCD on the isotropic lattices in the quenched approximation. 

5.5 Covariance of correlators 

Let us consider the pseudoscalar correlator: 

C PS (p,t) -- 



(O PS (x)O PS (0)) 
Z\p) 



2E(p) 



exp(-E(p)t) 



(large t) 



(5.7) 



with Z(p) = (0\O(x)\PS(p)). Here we employ the covariant normalization. For the local pseu- 
doscalar density operator 0(x) = q(x)j5q(x), if the Lorentz covariance is sufficiently restored, 
Z(p) does not depend on the momentum p. Then 



R(p) 



E(p)Z{pf 
m PS Z(0) 2 



(5- 



probes the restoration of covariance as the deviation from unity. In figure 12, we show the 



momentum dependence of R(p) measured on j3 = 5.95 and 6.10. At each j3, the quark mass is 
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Figure 12: The covariance of the pseudo-scalar correlators. The left figure is at /3 = 5.95 and 
k = 0.1230, and the right one is at (3 = 6.10 and n = 0.1230. 

the largest one used in the light hadron spectroscopy. We find that R(p) at finite momentum 
is consistent with R(p = 0), while higher momentum states suffer from large statistical fluc- 
tuation. This feature is particularly important in the calculation of form factors, in which the 
finite momentum states play an essential role. 



5.6 Systematic errors of the spectrum from calibration 

To estimate the systematic effect due to the uncertainty of calibration, we obtain the spectrum 
at the same k's with slightly shifted bare anisotropy, j' F = 7^ + 5"fF- We set 5"fF = 0.1, which 
implies about 2.5 % shift of bare anisotropy. Figure |l^ shows the result for shifted jf together 
with the result for j F , for (3 = 5.75 and 6.10. There are small systematic downward shifts 



in the fitted lines. The spectrum at the physical quark masses are listed in Table 13 as the 
dimensionless combination, m^rQ and m/mx*- The difference between the masses with 7^ 
and 7JJ1 is slightly amplified toward the chiral limit. Even for the lightest mass in each species, 
the difference is at most around 1 %. This implies that the uncertainty of hadron masses at 
the physical (u,d) and s quark mass are about half of uncertainty in 7^. With the relativistic 
dispersion relation, 7^ at m q = has been determined at each (5 within about 2 % ambiguity: 
the statistical error of 1 % and the systematic error of 1 % in the form of fit. Therefore there 
is 1 % level uncertainty in the hadron spectrum due to the uncertainty in calibration. This 
feature would make the anisotropic lattice promising for future physical applications. 



6 Conclusion 

In this paper, we studied the 0(a) improved quark action on the anisotropic lattice with 
anisotropy £ = a a /a T = 4. The bare anisotropy 7^, with which £p = £ holds, is determined 
for the whole quark mass region below the charm quark mass, including the chiral limit in 1 
% statistical accuracy. In the massless limit, there is also about 1 % systematic uncertainty in 
extrapolating 7^ to m q = 0. 

The uncertainties in the calibration due to the discretization errors are studied by changing 
the physical inputs or conditions, (i) We have shown that the dispersion relations for the 
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Figure 13: The spectra with shifted jf (filled symbols) together with the results at 7^ (open 
symbols). The solid lines and dashed lines represent the fit results at 7^ = 7^ and 7^ + ^7^, 
respectively. 



Table 13: Spectra with shifted 7^ in the dimensionless combinations, 
parameter J is also listed. 



The dimensionless 







m^rQ 






m/m* K 







5.75 


5.95 


6.10 


5.75 


5.95 


6.10 


p 


2.084(26) 


2.121(40) 


2.239(28) 


0.891(11) 


0.890(17) 


0.898(11) 


K* 


2.322(21) 


2.357(32) 


2.448(24) 


- 


- 


- 


<f> 


2.560(17) 


2.592(24) 


2.658(19) 


1.109(12) 


1.110(17) 


1.102(12) 


N 


2.946(37) 


2.913(52) 


3.019(43) 


1.260(16) 


1.223(22) 


1.212(17) 


A 


3.166(32) 


3.161(45) 


3.236(38) 


1.363(15) 


1.341(21) 


1.323(16) 


£ 


3.216(33) 


3.199(46) 


3.282(39) 


1.383(15) 


1.356(22) 


1.340(17) 


77 


3.486(29) 


3.485(41) 


3.545(35) 


1.507(15) 


1.490(22) 


1.468(17) 


A 


3.514(66) 


3.614(96) 


3.808(96) 


1.503(28) 


1.517(41) 


1.527(38) 


S* 


3.748(58) 


3.847(85) 


4.014(85) 


1.610(26) 


1.625(37) 


1.627(34) 


7T* 


3.983(50) 


4.080(73) 


4.221(75) 


1.717(24) 


1.734(34) 


1.728(30) 


n 


4.217(43) 


4.313(62) 


4.427(65) 


1.824(23) 


1.842(32) 


1.829(27) 


j 








0.3838(43) 


0.3871(86) 


0.3614(42) 



pseudoscalar and the vector mesons give values of 7^ which differ by 1% at (3 = 5.75, while 
they show no difference at (3 = 5.95 and 6.10. (ii) Two different choices of the lattice dispersion 
relations, namely the naive continuum form and the Klein-Gordon form, also lead to the results 
which differ by 3% for (3 = 5.75, but we found no difference at (3 = 6.1 with m q < 0.2a" 1 . 
These systematic uncertainties tend to vanish toward the continuum limit. 

The light hadron spectrum was studied using the central value of tuned bare anisotropy, 



tM 



m n 



0). We found that it is consistent with the result on the isotropic lattice by UKQCD 



Collaboration. It was found that a change of 7^ by 2 % would lead to a change of the spectrum 
by 1 % for the physical quark masses. We also investigated the Lorentz invariance of the matrix 
element of the pseudoscalar operator as a consistency check. 
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The main disadvantage in using the anisotropic lattice would lie in the additional systematic 
uncertainty caused by the calibration. There are two types of errors in j F . The first type 
consists of the statistical error and the error in the chiral extrapolation, which was estimated 
to be 2 % level. The second type consists of O(aa) and 0(a 2 ) systematic uncertainties, which 
was estimated to be 4% at /3 = 5.75 and smaller for larger j3. In total, there is 6 % ambiguity 
at j3 = 5.75, which corresponds to our coarsest lattice, and 2-3 % ambiguities at (3 = 6.10. 
The relative errors in hadron spectrum are half of those in j F . Since the contribution from the 
second types of error vanishes in the continuum limit, we expect to obtain the hadron spectrum 
to 1 % accuracy in the continuum limit. This result is encouraging for further applications. 
The anisotropic lattice would already be applicable to quantitative studies which requires a few 
percent accuracy. To achieve higher accuracy, nonperturbative tuning of the clover coefficients 
are required. 

Since the range of quark masses where the systematic errors are under control covers the 
charm quark region, it is also important to apply the present anisotropic lattice simulation to 
the charmonium and D meson systems. 
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